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Abstract. The SCUBA submillimetre camera has opened up new pos- 
sibilities for tracing the evolution of active star formation in dusty galaxies 
^N| ■ to high redshift, with profound implications for our understanding of the 

K*" \ star formation history of the Universe and the contribution from galaxies 

^-p ■ to the anisotropy of the microwave sky. We review results from several 

submillimetre surveys started during SCUBA's first year of operation, 
ps) ■ and discuss their future development, together with other projects that 

<^ \ will greatly improve our understanding of the extragalactic point source 

0^ ■ contribution to the submillimetre sky in the era of MAP and Planck . 

O ' !• Introduction 

^ ■ Advances in modern observational astronomy are usually led by innovation in 

detector technology, rather than through inspiration on the part of observers. A 
good example of this is the recent upsurge of interest in submillimetre astron- 
KA , omy, following the installation of the Submillimetre Common User Bolometer 

'i_j ■ Array (SCUBA0: Holland et al. 1999) on the James Clerk Maxweh Telescope 

on Mauna Kea, Hawaii. The (sub)millimetre is one of the few regions of the 
electromagnetic spectrum remaining largely unexplored by modern astronomy, 
and, as we detail below, SCUBA is the first instrument to cross the sensitivity 
threshold beyond which serious extragalactic astronomy in the submillimetre 
becomes not only possible, but very powerful. This has led several groups to 
start survey programmes with SCUBA, and this article reviews their progress 
to date. These surveys impact on the study of CBR foregrounds in three ways: 

1. the central wavelength of the most sensitive SCUBA band (850/im) matches 
that of one of the channels of the Planck High Frequency Instrument (HFI: 
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Puget et al. 1998), so SCUBA is directly sampling one of the foregrounds 
to be seen by Planck, albeit at fainter flux densities than HFI will reach 

2. the source counts resulting from the SCUBA surveys will provide impor- 
tant constraints on the models of galaxy evolution required to predict the 
level of point source contamination in microwave sky maps produced by 
CBR experiments in less directly studied passbands 

3. SCUBA observations may resolve the cosmic infrared/submm background 
(Puget et al. 1996) detected in COBE maps. 

These issues are addressed more directly in the contributions to this volume 
by Eric Gawiser, Bruno Guiderdoni, Luigi Toffolatti, and collaborators, so we 
shall concentrate here on observational aspects of SCUBA survey programmes. 

2. Why study galaxies in the submillimetre? 

IRAS revealed a strongly evolving population of starburst galaxies forming stars 
at rates of a few hundred solar masses per year in heavily obscured regions, 
reminiscent of the formation of massive stars in giant molecular clouds in the 
Galaxy. Typical starbursts emit ~ 50% of their bolometric luminosity in the 
rest-frame far-infrared, as the UV radiation from massive young stars heats the 
dust around them, while this fraction can increase to -^ 90% in ultraluminous 
infrared galaxies (ULIRGs: those with LpiR > IO^'^Lq), although some of these 
are partially powered by active galactic nuclei (Genzel et al. 1998). The absorp- 
tion of such a large fraction of the UV/optical radiation from massive stars in 
starburst galaxies necessarily means that studies performed in those bands will 
underestimate the star formation rates in individual starbursts (as indicated by 
Cram et al. 1998, from the comparison of estimates of the rates of star formation 
in a sample of ~ 700 local galaxies derived from U band, Ha, far~infared and 
decimetric radio continuum luminosity^, as well as allowing the possibility of 
the complete omission of a population of heavily-obscured starbursts from the 
cosmic star formation census. Since most massive stars are made in starburst 
galaxies this can have a significant effect on the global star formation rate, and, 
once extinction corrections (Meurer et al. 1997, Pettini et al. 1998) are applied, 
the star formation history of the Universe deduced from UV/optical studies can 
change substantially (e.g. Steidel et al. 1999) from the well-known shape of the 
much-touted "Madau diagram" (Lilly et al. 1996, Madau et al. 1996). 

In Figure 1 we plot the spectral energy distribution (SED) of the prototyp- 
ical local starburst galaxy M82, which shows the characteristic starburst shape 
of a peak at ~ 60 fim, and a steep decline to longer wavelengths, well fitted by 
a single temperature greybody curve with emissivity oc zv^~^ and temperature 



■^While these and most other conventional methods of estimating the star formation rate of 
a galaxy determine the luminosity from the massive stars it contains, the exact stellar mass 
range probed by each method is slightly different, and the conversions from luminosity density 
to star formation density required to put each point on the "Madau diagram" are uncertain 
to factors of a few (e.g. Rowan-Robinson et al. 1997) due to differing assumptions about the 
stellar initial mass function. 
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Figure 1. The SED of the prototypical starburst galaxy M82 (from 

Efstathiou et al. 1999: ERRS). Observational data points are plotted, 

together with model SEDs resulting from two bursts of star formation 

(dotted and dashed lines), whose sum (solid line) provides a good fit 

to the galaxy's SED over the range 0.5-850 /xm: see ERRS for details. 
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Figure 2. Starburst submillimetre flux density as a function of z. 
The solid (dotted) line shows the result of taking M82 (the solid line 
of Fig. 1) and observing it at 850//m at increasingly high redshift in a 
Universe with f^o = 1 (^o = 0.1) and A = 0. The dashed (dot-dashed) 
lines show the results of doing this at 350/um (2mm), while the dot- 
dot-dot-dashed line shows the effect on the solid line of replacing M82 
(L ~ IO^'^Lq) with Arp 220, which is about ~ 60 times more luminous, 
and has a slightly different SED (A. Efstathiou, priv. comm.). 



~ 50 K (although the physics of the dust emission is certain to be more comph- 
cated than this simple model might seem to imply). The steepness of this decline 
means that such galaxies have high negative fc-corrections in the submillimetre 
- i.e. as they are observed in a fixed passband at increasingly high redshift, the 
rest-frame wavelength of the observed emission moves towards the peak and so 
the rest-frame flux increases - which counterbalance the inverse square law effect 
due to the increasing distance, with the result that the observed submillimetre 
flux (in a fixed passband) of a starburst galaxy of a given luminosity is almost 
independent of redshift for 1 < z < 10 or so (for JIq = 1 - the effect depends on 
cosmology through the luminosity distance), as shown in Figure 2. 

3. The Submillimetre Common User Bolometer Array (SCUBA) 

SCUBA consists (Holland et al. 1998) of two arrays of bolometers (37 in the long 
wave array, 91 in the short wave array) allowing simultaneous imaging in pairs 
of submillimetre bands (the most sensitive being 850;Um and 450/Um) over a field 
of view of ~ 2.3 arcmin diameter, or photometry at 1.1, 1.3 or 2.0 mm using 
single detectors in each band. The bolometer feedhorns are close-packed, but 
finite gaps between them necessitate a series of offset observations (performed by 
"jiggling" the secondary mirror of the JCMT) to sample the field of view fully. A 
total of 64 pointings is required to produce fully sampled maps at both 850 and 
450 ;um, and these have diffraction-limited angular resolutions of 14.7 and 7.5 
arcseconds (FWHM), respectively. Each bolometer is ~ 100 times more sensitive 
than the previous JCMT single-element broadband photometer (UKT14) which, 
combined with its large number of pixels, means that SCUBA reaches the plateau 
in the flux-redshift plot for luminous starbursts in a sufficiently short time and 
large field of view to make large submillimetre surveys feasible for the first time. 

4. Submillimetre surveys with SCUBA 

Four groups have reported results from SCUBA surveys, and several others are 
underway. Two complementary strategies have been adopted: (i) simple map- 
ping of blank fields; and (ii) pointed observations of massive clusters of galaxies, 
using their gravitational lensing effect to probe deeper than otherwise possible in 
a given integration time. Harnassing gravitational lensing amplification probes 
deeper (simplifying follow-up observations in other wavebands, as well as initial 
source detection with SCUBA) but it can introduce into quantitative results 
uncertainty resulting from imprecision in the cluster lens model, and, given the 
large beam size, there is the possibility of misidentifying emission from cold dust 
around the central galaxy in the cluster core - as seen by IRAS in some clusters 
(Bregman, McNamara & O'Connell 1990) and as expected in cooling fiow mod- 
els (see Fabian 1994 for a review) - with emission from a gravitationally-lensed 
background galaxy along the line of sight. 

4.1. The SCUBA Lens Survey 

The first SCUBA survey results were reported by Small, Ivison & Blain (1997) 
who detected six sources in observations of the two clusters A370 and CI 2244-02 



at z = 0.37 and 0.33 respectively. They deduced a source density of (2.4 ± 1.0) x 
10^ deg~^ at a 50% completeness limit of ~ 4mJy at 850 ;um, requiring strong 
evolution of the starburst population at z > 1 (on the basis of the IRAS 60 /um 
luminosity function of Saunders et al. 1990, and the galaxy SED models of Blain 
& Longair 1996). The SCUBA Lens survey is now complete (Small et al. 1999), 
with 17 sources detected at 3a or better, from observations of seven clusters 
over the redshift range 0.19 < z < 0.41, covering a total area of 0.01 deg^ to a 
depth of 2mJy. Notable amongst the follow-up of the sources is the study of the 
galaxy SMM J02399-0136 (Ivison et al. 1999) detected in the field of A370 and 
identified with a z = 2.8 galaxy. It exhibits a hybrid of features associated with 
star formation (high molecular gas mass deduced from CO emission, high radio 
and Ha luminosities) at a rate of several thousand solar masses per year, and 
with the presence of an active nucleus (optical emission line properties and radio 
morphology): this mix of features is also seen in some hyperluminous infrared 
galaxies (e.g. F10214+4724: Lawrence et al. 1994). 

4.2. The UK Submillimetre Survey 

Our own UK Submillimetre Survey comprises two separate blank field pro- 
grammes - a narrow, deep survey and a wide, shallow survey ~ which will 
combine to produce a dataset capable of tracing the form of the evolution of 
the starburst population to high redshift. We have completed our narrow, deep 
survey (Hughes et al. 1998) which was undertaken in the Hubble Deep Field 
(HDF: Williams et al. 1996) during a period of exceptionally good weather con- 
ditions. We covered an area of ~ 5.5 arcmin^ to a la noise level of 0.45 mJy 
at 850 /xm, putting us beyond the classical confusion limit, given the size of the 
SCUBA beam. This complicated the identification of discrete sources, but, on 
the basis of simulations (see Hughes et al. 1998 for details), we conservatively 
accepted as real and discrete a sample of 5 sources with fluxes above 2 mJy. 

The principal reason for undertaking our narrow, deep survey in the HDF 
was the existence of a wealth of data in other wavebands in that field, which 
would help with the process of identification of our SCUBA sources. For our five 
sources we found four plausible associations, with galaxies at redshifts (spectro- 
scopic or photometric) z ~ 1, 2, 3 and 3, while the fifth source had two candidate 
counterparts, one with a photometric redshift of z ~ 4 and the second a faint 
galaxy lying in the Hubble Flanking Fields and therefore not amenable to photo- 
metric redshift estimation. Note that the optical/far-infrared flux ratios typical 
of starburst galaxies do not preclude our SCUBA sources being high-redshift 
galaxies too faint to be detected in the optical HDF, despite its unprecedented 
depth, and also that only our brightest source (with a flux of 7 mJy at 850 
^m) is detectable at 1.3 mm with the IRAM interferometer, which would yield 
a more accurate source position, simplifying its identification. 

4.3. The Hawaii survey 

In the same issue of Nature in which we reported on our HDF survey, Barger et 
al. (1998) presented results of pointings in two blank fields: one in the Lockman 
Hole, and one in Hawaii deep survey field SSA13. They detected two sources 
above 3 mJy, the brightest of which they deduce to lie in the range 1.5 < z < 3.5. 




Figure 3. Integral source counts from the four 850;uni surveys de- 
scribed in the text: empty circles (Small et al. 1999); filled circle 
(Small et al. 1997); empty triangle (Barger et al. 1998); filled triangle 
(Eales et al. 1999); and filled squares (Hughes et al. 1998). For clarity, 
filled symbols are plotted O.lmJy brighter than their true flux limit: 
the line - iV(> S) = W'^{S/lm3y)-^/'^[deg-^] - is not a fit to the data. 



4.4. The Canada UK Deep Submillimetre Survey 

A second wide-area survey is being undertaken in two of the fields of the 
Canada-France Redshift Survey (CFRS), which total 200 arcmin^. Results from 
the first 22 arcmin^ have been presented by Eales et al. (1999), who report de- 
tections of 12 sources with S'sso^m > 2.8 mJy. Lilly et al. (1999) found within 
the CFRS associations for six of these sources, while a further two have tenta- 
tive identifications to the CFRS photometric limit of /ab = 25: of these, four 
have spectroscopic or photometric redshifts below z = 1, four lie in the interval 
1 < z < 3, while the four unidentifed sources are assumed to have z > 3. 



5. Summary 

The results reported to date appear to present a fairly consistent picture. Fig- 
ure 3 shows that the integral source count distribution is already becoming con- 
strained over more than an order of magnitude in flux density, while Lilly et al. 
(1999) note that the redshift distribution they derive appears broadly consistent 
(given the small numbers of sources, and the uncertainty in estimating photo- 
metric redshifts at z > 1) with that of our HDF sources, and that the /ab band 
magnitude distributions of the galaxies associated with the CFRS and SCUBA 
Lens Survey sources also appear to agree well. All four groups broadly ascribe 
their results to a population of luminous starburst galaxies, which resolve most 
(perhaps nearly all?) of the cosmic far-infrared/submillimetre background, have 
SEDs broadly similar to that of Arp 220, and are observed to be forming stars 
at rates of several hundred solar masses per year. Much of this activity is taking 



place at redshifts where previous UV/optical studies (Madau et al. 1996) sug- 
gested there should be little star formation, although uncertainty over the local 
template to use for these galaxies means (Eales et al. 1999) that it is not clear 
whether these two views of the star formation history really are inconsistent. 
The revision upwards of the UV/optical star formation rates at z > 2 (Steidel et 
al. 1999) does, however, make a consistent picture emerge, in which luminous, 
dusty starbursts at high z produce a large fraction of the Universe's stellar mass. 

It is tempting to identify these objects as the long-sought precursors of 
today's luminous elliptical galaxies, and perhaps this will become clear by the 
time these surveys are complete. The two wide area surveys described here 
should detect a few hundred sources in total, sufficient to characterise the nature 
and evolution of this population well (although optical spectroscopy of most of 
these objects will not be possible, even with 10m class telescopes). In particular, 
our UK Submillimetre Survey has been awarded the time needed to complete 
our shallow {Sa depth of 8 mJy at 850^m) survey of 640 arcmin^ in three fields 
with deep ISOPHOT 175 //m data: these sources should be bright enough to 
follow up at 1.3mm with the IRAM interferometer, to produce more precise 
positions, making the association procedure easier, while non-detection of any 
sources in the deep 175/im data robustly indicate that the galaxy lies at z > 2. 

As discussed in other contributions to this volume, the SCUBA results de- 
scribed here are already providing information on the extragalactic point source 
population that will be important to MAP, while, by the time Planck is launched 
in 2007, there will be several other very powerful facilities available for studying 
galaxies in the (sub)millimetre. Firstly, SCUBA itself may be upgraded^, with 
suggested improvements capable of increasing sensitivity so that the confusion 
limit at 850/im can be reached in 2-3 hours of good weather, rather than ~20 
as at present. The first of the completely new facilities is the Large Millimeter 
Telescope (LMT^, a 50m single dish telescope to be built by a US/Mexican 
consortium at a site in Mexico at an altitude of 4,600m. The LMT will operate 
over the range 0.85-3. 4mm, and its Bolocam broad band camera should reach 
a lOo" sensitivity of 0.3 mJy in 1 hour at 1mm (D. Hughes, priv. comm.). Most 
important is the planned Large Millimetre Array (LMA), which is likely to re- 
sult from a merger of ESO's Large Southern Array|^ project and the Millimeter 
ArrayQ planned by NRAO for the same exceptional site at Llano de Chajnantor 
in Chile. The design of the LMA is still under discussion, but it will provide the 
possibility of sub-arcsec resolution (sub)millimetre imaging by interferometry to 
very faint levels, albeit in a small field of view. This would allow detailed follow 
up of sources found in surveys made with SCUBA, LMT, FIRST^, etc, and also 
with the HFI on Planck itself. It is clear that the continued rapid development 
of (sub)millimetre astronomy over the next decade will provide the greatly im- 



see www . j ach . hawaii . edu/ JCMT/scuba/upgrades/up . html 
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proved understanding of the extragalactic point source population required for 
the success of future CBR experiments hke MAP and Planck as weh as facihtat- 
ing important secondary science with these missions, studying the astrophysics 
of the point sources they wih remove from their maps of the microwave sky. 
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